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Abstract 

A comparative study of the liquid-drop model (LDM) type expansions of energy 
E and compression modulus Ka is made within the energy density formalism using 
Skyrme interactions. As compared to the energy expansion, it is found that, in the 
pure bulk mode of density vibration, the LDM expansion of Ka shows an anomalous 
convergence behaviour due to pair effect. A least squares fit analysis is done to 
estimate the minimum error, one would expect even with synthetic data due to 
the inherent nature of the LDM expansion of Ka as well as the narrow range of 
accessible mass number A, in the values of the various coefficients. The dependence 
of the higher-order coefficients like curvature and Gauss curvature on the coupling 
f3 c between the bulk and surface parts of the monopole vibrations is analytically 
studied. It is shown that the Ka— expansion including the dynamical effect ( A— 
dependence of (3 C ) shows an 'up-turn' behaviour below mass number about 120, 
suggesting the inapplicability of the LDM expansion of Ka over this mass region. 



PACS numbers: 21.65.+f, 21.10.-k 



I. INTRODUCTION 



Nuclear matter incompressibility Koo is of fundamental significance as it is an impor- 
tant ingredient in the nuclear equation of state, which influences several astrophysical 
phenomena as well as high energy heavy-ion collision processes. The nuclear breath- 
ing mode or isoscalar giant monopole resonance (GMR) is the most promising source of 
information on K^. 

However, the nuclear breathing mode actually determines the finite nuclear compres- 
sion modulus Ka-i and not K^. Hence, to determine the value of K^ one has to extrap- 
olate from finite nuclei to nuclear matter, which is a highly non-trivial task. Inspired by 
the success of the nuclear mass formulae, a liquid-drop model (LDM) like approach is 
currently being adopted to extract Koo from its finite nuclear value Ka- In this approach, 
one first expresses the finite nuclear compression modulus Ka in terms of the experimental 
breathing-mode energies E gmr using the relation, 



where the experimental values of the mean square radius < r 2 > are used; and then one 
supposes an LDM expansion of Ka to be valid, in analogy to the semi-empirical mass 
formula; 



where (3 is the asymmetry parameter. The volume coefficient K v determined from a fit to 
the available breathing mode energies is identified with K^, which is strictly valid only 
in the scaling model. 

Following this approach, Sharma et al Jl|, 0] determined K v to be about 300 MeV 
using their data on various Sn and Sm isotopes, and on those of 24 Mg, 90 Zr and 208 Pb, 
which was later contested by Pearson || . Subsequently, Shlomo and Youngblood HQ made 
an extensive least-squares fit analysis by taking into account the available set of data on 
breathing-mode numbering about 46, spread over a mass region 28 < A < 232. They 
concluded that this complete data set is not adequate to determine the value of Koo to 



K A = m < r 2 > E 2 mr /h 2 




K A = K V + K S A~^ 3 + K C A~ 2 / 3 + K Cou Z 2 A-^ + K p (3 2 



(2) 
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better than a factor of about 1.7 [ 200 -350 MeV ]. Further, they also observed in their 
analysis that a free least-squares fit to Eq.(2) leads to errors exceeding 100% in all the 
coefficients. As observed by Shlomo and Youngblood, such large errors may arise partly 
due to the break-down of the scaling approximation over light nuclei, deformation effects 
etc. Another reason may be that the presently available data are inadequate for fixing 
all the parameters in Eq.(2). We feel that, possible slow convergence and correlations 
among the coefficients, will also contribute toward this error. However, it is not possible 
to pinpoint what fraction of the error arises out of the quality and inadequacy of the data, 
and that which arises due to the nature of convergence of the LDM expansion of It 
is therefore desirable to make an extensive theoretical study of the nature of convergence 
and the effect of correlations to arrive at a conclusion. 

Studies in this regard have been attempted [5-9] over the years and have normally 
been done in the scaling model, since an analytical relation between the experimental 
breathing-mode energies and Ka is available only within this model. But, it is known 
from the works of Brack and his collaborators[10-12] that there can exist several types of 
density vibration, the scaling mode being a particular one. It was also shown by them 
in a hydrodynamical approach that, the degree of coupling between the bulk and surface 
parts of the density vibrations, sensitively depends upon the mass number A. 

Let us consider the ground-state density of a nucleus to be given by a Fermi-function 
p(r) = p /[l + e^~^] where p a , R and a are the central density, half density radius and 
diffuseness paramters respectively. When it undergoes monopole vibrations, the density 
function will be compressed or decompressed. In this state, it is assumed that the den- 
sity function can still be represented by a Fermi-function. Then, following Brack and 



Stocker pJJI , one has 

P c (r) = 

1 + e 



P c (r) = \_ R) (3) 



where p c and a c are defined as 



Pc = PoQ 

a c = aq p ° = a(p c / Po f c (4) 



The paramter q represents the amount of compression or decompression, j3 c gives the 
degree of coupling between the bulk and surface parts of the density vibration and a c 
is the corresponding value of the surface diffuseness. The scaling property of the half- 
density radius R is defined by the number equation given by Eqs.(6,12). In literature [10- 
12], the mode of vibration pertaining to f3 c = —1/3 is normally referred to as scaling 
mode of vibration. In realistic hydrodynamical calculations [|Tl~| , [12"|1 , which describe the 
experimental data on breathing-mode reasonably well, it is found that (3 C can vary from 
-0.23, in the case of 208 Pb, to -0.84 for 90 Zr. Further, f} c tends to ±oo for lighter nuclei 
like 48 Ca \j3 c = -24] and 40 Ca [f3 c = +9.5], i.e. pure surface like mode. 

More importantly, when one makes use of an LDM expansion of Ka to fit the experi- 
mental GMR data, the various finite-size coefficients in the expansion(2), such as surface 
K s and curvature K c , must be A— independent, as in the case of standard nuclear mass 
formulae. Due to the mass number dependence of the dynamical coupling paramter /3 C , 
there is still a residual mass dependence in these coefficients obtained theoretically. The 
asymptotic value of K s obtained in the limit A — > oo is about —2K V , and corresponds 
to the pure bulk mode /3 C =0. This large, negative value of K s shall lead to a slower 
converging Ka— expansion as compared to the one obtained in the scaling model, where 
K a ~ —K v . Therefore, it is of interest and relevance to study the convergence properties 



of the LDM expansion of Ka for this particular mode. This has been stated|jTTJ to be a 
well converging series, which may not be conclusive as we shall see. 

In view of the above discussions, we make an analytical study of the LDM expansion 
of Ka derived with /3 C =0. It is found that this particular expansion shows an anomalous 
convergence behaviour in the sense that certain higher-order terms are equal in magni- 
tude and opposite in sign to the preceding lower-order terms in the expansion, resulting 
in, what is termed here as, 'pair effect'. Secondly as the LDM expansion of Ka is rela- 
tively new compared to the well established energy expansion, it is worthwhile to make 
a comparative study of both these expansions at a fundamental level. Thirdly, we also 
make a least-squares fit analysis to get an idea regarding the minimum error in the val- 
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ues of the various coefficients, one would expect even with synthetic data because of the 
inherent nature of the LDM expansion as well as, the inadequacy of the narrow range 
of A to fix the higher-order coefficients. Finally, to ascertain the goodness of the Ka— 
expansion for the extraction of the various coefficients, we have studied the the f3 c depen- 
dence of curvature K c and Gauss curvature K Q coefficients, and thereby extending the 
pocket model calculations to higher-order coefficients. It is found that, the validity of the 
LDM expansion of Ka in particular cases such as scaling and pure bulk mode, does not 
necessarily quarantee unambiguous extraction of K^ in realistic situations. This is shown 
by taking the dynamical effect (the A— dependence of f3 c ) into account, where the nuclear 
incompressibility Ka is found to exhibit an 'up-turn' behaviour below mass number about 
120, which may be of non-leptodermous origin. 

II. THE MODEL 

Here, we analytically derive the LDM expansion of nuclear incompressibility Ka start- 
ing from first principles, and using the concept of leptodermous expansion of energy. The 



model considered here is essentially the same as the one given in Ref.|13|, which we have 
now generalised to consider various modes of monopole vibrations. 

We consider a symmetric system of A fermions, without the Coulomb interaction, 



whose total energy E, in the framework of energy density formalism ||14|| is given by, 

E = J e [p(r)]d 3 r (5) 
with the number constraint equation, 

A = f p(r)d 3 r (6) 



where, 

e(r) = —r(r) + v(r) (7) 



For the kinetic part r(r), we use the extended Thomas- Fermi functional [15| upto 
0(h 2 ), which for a symmetric system is 



5 00 p 6 



For the potential part v(r), we use the standard Skyrme forces, without the spin-orbit 
contribution, 

v(r) = — p + -p + ~^^rp + — g^(Vp) (9) 

Since our objective is to make an analytical study, we do not attempt to solve the 
Euler-Lagrangian equations for self-consistent densities. Further, to systematically study 
the general nature of the LDM compressibility coefficients and their dependence on the 
coupling paramter /3 C , it is desirable to parametrise the density function p(r) to be a Fermi 
function. Therefore, the ground-state density function p(r) for a nucleus is given as 

P(r) = P \ r . R) (10) 
1 + e « 

where p a , R, and a are the central density, half-density radius and the diffuseness pa- 
rameter respectively. In order to obtain analytical expressions for the various finite-size 
compressibility coefficients in terms of the coupling paramter /3 C , we define the general 
density distribution corresponding to a compressed/decompressed state to be of the form 
given by Eq.(3). Thus, our whole derivation is of general nature, and the ground-state 
properties can be determined by imposing q = p c j p a = 1. 

Normally, one expands the finite nuclear ground-state energy into volume, surface etc 
using Taylor's series. Then, each coefficient in the expansion is calculated(for e.g., see ||) 
making use of the semi-infinite nuclear matter approximation . Instead, we here arrive 



at the LDM expansion of Ka using the generalised Sommerfeld lemma fl6 , 17] . One can 
then systematically express the total energy as a sum of volume, surface, curvature and 
Gauss curvature. Further, as we are interested in the convergence properties of the LDM 
expansion of Ka, we go upto Gauss curvature order in both the energy and compressibility 
expansions. 

Then, the total energy E at any state of compression/decompression given by Eq.(5) 
results in the leptodermous expansion of E as, 

4irR 3 

E = E V + E s 4nR 2 + E c 8ttR + E 4n (11) 



The terms E v , E s , E c and E Q are the volume, surface, curvature and Gauss curvature 
contributions to the total energy E. Then, the half-density radius R can be given in 
terms of p c , by using Eq.(6) as 



R ~ r nr A 1/3 



3rlA 2 / 3 



(12) 



where -^j^Pc = 1- Subtituting the above expression for R in Eq.(ll), and retaining terms 
upto the order of A^ 1 , we obtain 

E 



el + eM- 1/3 + e*A^ 3 + eZA' 1 



(13) 



The expressions for the different coefficients in the above equation are 



e* v = PcPc 2/3 + PoPc + Pi«p c 5/3 + VzPc a+l 



(14) 



= -(367r) 1 / 3 a c (0.759p aPc - ^p c 1/3 + p p c 4/3 + 1.359a Pl p c 2 



6a 2 c 



(15) 



e„ = 



4a c 2 (6 7 r 2 ) 1 /3((i. 5 17 - ^) Pa p^ + ||p 2 / 3 + (1.973 



7T 



)pi«p 



,7/3 



7T 



+(0-5<& - y)p 3 p 



.+5/3) 



(16) 



-47ra c 3 ((2.602 - ^0.759)p a ^/ 3 + ^*§p c + (4.423 - ^1.359) Pl a Pc 8 / 3 
3 b ai 3 



7T" 

"3 



2tt_ 2 
3 



(17) 



where p 

3/ 



n 2 



3to 



*3 ^ — 3ti+5fa „ 
16 > Pi - 16 ' ^2 



and a 



5 -(1.57r 2 ) 2 / 3 , P — and 7 = §• The quantities gg s , <7g c and are calculated using the 
integral 



k I „,k 



lo 



1 + -ire 



k„—uv 



(1 + e-") 1 ' 

where v = a + 2, and = 77 (°\ <& = 2?^ and <£ = 
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It must be noted that all the coefficients in Eq.(13) are function of the ground-state 
values of central density p Q and the diffuseness parameter a, and also are function of the 
compression variable q and the coupling paramter f3 c . p Q and a are normally determined 
using the energy minimisation criteria. Due to the finite compressibility of the nuclear 
fluid, both the quantities could be functions of the size of the system, characterised by 
their mass number A. Therefore, Eq.(13) can not be identified as the LDM expansion of 
the nuclear energies. 

To extract the A-dependence of the energy coefficients in Eq.(13), we define p {A) = 
Poo + 8p, where p^ is the symmetric infinite nuclear matter saturation density. p (A) can 
be determined using the saturation condition, 

d(E/A) 



= (19) 



dp 

Remembering that the above derivative is a total one, we define 

ME /A) ME/A) 
Po I-^ll^ = [^Hl*- ( 2 °) 

The saturation condition given by Eq.(19) now becomes, 

e' v + e' s A- 1 '* + e' c A- 2 l* + e^ 1 + e'^ 4 / 3 + e^ 5 / 3 = (21) 

Here, the prime denotes total differentiation of e*'s with respect to the central density 
variable. To determine all the coefficients at the density of symmetric infinite nuclear 
matter , we use the definition 

Po = Poo + Sp (22) 

where 5p contains all the possible finite-size effects. 

Making Taylor's expansion of all the coefficients in Eq.(21) around upto 0(5p) 2 , 
we obtain, 

eU" 1/3 + e'A~ 2 / 3 + e'A- 1 



5p = 



< + e(M-V3 + e£M- 2 /3 + e'^A' 1 + 0.5(e;" + efA' 1 / 3 + • • -)Sp 



(23) 
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As the above equation contains 5p on the right hand side also, we determined the same 
by an iterative method. The zeroth order expression of 5p is equal to the above equation 
without the bracketed term in the denominator. 

To arrive at the LDM expansion of E/A and K^, we should be able to expand 5p in 
powers of A -1 / 3 . Unlike in Ref. [13|, here we consider the influence of higher-order terms 
like e", e" etc by making a binomial expansion of the denominator in Eq.(23). Considering 
terms upto third order, we have, 

1 + (fl _ e " & s W-l/3 + , , 

~ 1 + f\(n)A- l l* + f 2 (n)A- 2 / 3 + ■■■ (24) 



where, 



fi(n) 



—n 



—n 



—n 



pin P i 
2«) 2 

2«') 2 



e"'e' 

ll\2 



2 e: 



piiipi'pi 
2?ZP 



//// / 2 

6«) 3 



+ 



nfn + r 



e'V 
2«) 2 



It may be recalled here that prime denotes total differentiation with respect to the cen- 
tral density variable for a fixed value of f3 c . One then has for any arbitrary function 
/ [p c ,a c (p c )], 

Po \dq) 
1 (£f 



f" = 



'dn 

dpcj 
'dPf 



\p - 



\q=\ 



dpi 



Pi 



1 9=1 



(25) 



A. LDM expansion of energy 

Now, the LDM expansion of E/A can be obtained in a straight-forward manner by 
performing Taylor's expansion of each of the term in Eq.(13) around as, 

| = e* v ( Poo ) + e* s { Poo )A- l l* + ■ ■ ■ + e'MA- 1 
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+ (e s \ Poo )A-^ + ■■■ + eJip^A- 1 ) bp 

+1 (e v "( Poo ) + e s "{ Poo )A-^ + ■■■ + e '\ Poo )A- 1 ) {5pf 

+1 (e v "\ P oo) + e s '"{ Poo )A-^ + ■■■ + e '"{ Poo )A-') {8pf + ■■■ (26) 

Then, by using Eqs. (23-24) in the equation (26), and grouping terms with same power of 
A, the complete LDM expansion of E/A upto 0(A~ l ) is 

E 

— = a v + a s A- 1/3 + a c A~ 2/3 + a Q A~ l (27) 

where the various LDM coefficients in the above equation are defined as, 

a v = e*(poo) 
a s = e*(poo) 

1 

a c = e*(poo) + (due' s ) + - (d 22 e") 
a = eo(poo) + (d n e' c + d 12 e' s ) 

+l(d 22 e' s ' + d 23 e^ + ±(d 33 e':) (28) 

It may be noted that unlike the leptodermous expansion given by equation (11), the 
LDM expansion of energy is an infinite series. This infinite nature comes due to the A^ 1 ^ 
expansion of 5p, which is derived by making a binomial series. Therefore, we feel it is 
essential to go upto atleast Gauss curvature order in the LDM expansions of energy and 
compressibility, while one studies their convergence properties. 
The explicit expressions for the various factors du,di 2 etc are; 

(i) d u = ~e' 8 

du = -^(e' c + eMl)) 

dis = -^(e'o + e'Ml)+e' s f 2 (l)) 

(ii) d 22 = —^e' s 2 

d23 = -^(2e' s e' c + e' s 2 h(2)) 
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(Hi) d 33 = -ije; 3 . (29) 

The important point to be noted regarding these factors is that they all have an explicit 
dependence upon K^, which in turn controls the convergence behaviour of the LDM 
expansion. However, in the case of energy, the two leading coefficients volume and surface 
are independent of these factors, and thereby K^. Only the higher-order terms like 
curvature and Gauss curvature are dependent on K^. 

B. LDM expansion of nuclear compression modulus 

In the following paragraphs, we derive the LDM expansion of K A . The finite nuclear 
compression modulus Ka is calculated using the definition, 



2l d 2 E/A\ d 2 E/A , 

k a = Vo ( -nr- J = 9-772- U (3°) 

\pc = Po 



dpi ). dq 



for a given value of (5 C . 

Making use of equations (27-28) in the above definition, we have 

K A = K* v ( Po ) + K:( Po )A-^ + K:( Po )A- 2/S + K^A- 1 (31) 

where, 

KKPo) = 9p^||l 

2„* 



dq 



U=i ; i = v,s,c, and o (32) 



Now, the LDM expansion K\ can be obtained in a straight- forward manner by performing 
Taylor's expansion of each of the term in Eq.(31) around p^ as, 

K A = K* v (p 00 ) + K:(p 00 )A~ 1 / 3 + --- + K* (p 00 )A- 1 

+ (k:\ Poo ) + x;'( Poo )A-V3 + . . . + k:\ Poo )a^) 5p 

+\ (K"(Poo) + Kl"{ Poo )A^ + ■■■ + K* "( Poo )A- 1 ) (5p) 2 

+1 (K:"'(Poo) + K'"{ Poo )A-^ + • • • + Kl"\ Poo )A- 1 ) (Spf + ■■■ (33) 
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Here, we have included all the terms upto 0(Sp) 3 , so that correct estimates upto Gauss 
curvature order i.e. 0(A^ 1 ), can be obtained, which is desirable in this study of the 
convergence behaviour of the LDM expansion of K A . 

Then, by using Eqs. (23-24) in the equation (33) and grouping terms with same power 
of A, the complete LDM expansion of K A upto 0(A~ l ) is 

K A = K V + K S A- 1/3 + K c A- 2/3 + K Q A- 1 (34) 

where the various LDM coefficients in the above equation are defined as, 

K v = K* v ( Poo ) 

k s = k;{ Poo ) + (dnK) 

K c = K* c (p 00 ) + (d ll K' s + d 12 K' v ) + ^(d 22 K > :) 
K D = K* (p 00 ) + (d 11 K' c + d 12 K' s + d 13 K' v ) 

+ l - {d 22 K" s + d^K) + \ (d 33 K':) (35) 

The explicit expressions for the various factors d n ,di 2 etc are as given in Eq.(29). It may 
be noted that unlike in the case of energy expansion, all the finite-size coefficients (K s , 
K c , K ) explicitly depend upon through the factors du, d\ 2 etc. 

III. RESULTS AND DISCUSSIONS 

In this section, we discuss the basic nature of an LDM expansion of energy and incom- 
pressibility. In addition, the dependence of the various compressibility LDM coefficients 
on the coupling parameter f3 c is also studied. 

A. Structure of the LDM coefficients and dependence on (3 C 

The first and the important point is that in the case of energy, the bulk part totally 
decouples itself from all the surface effects due to the saturation condition e^(poo) = 0, 
in the leading order Sp, which is evident from Eq.(26). This, in turn, results in the two 
important LDM energy coefficients a v and a s , given in Eq.(28), being pure in nature, 
in the sense that there are no contributions arising from 5 p. In contrast, in the case of 
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compressibility, the volume gets strongly coupled to all the surface effects even through 
the leading order in 5p, and hence, the surface coefficient K s has an additional term 
duK' v , which explicitly depends upon through di%. This fact contributes to the 
essential difference between the LDM expansion of E/A and As it can be seen from 
Eqs.(28) and (35), it also gives rise to additional terms in the higher-order coefficients in 
the expansion of Ka- With this observation, we now calculate the values of the various 
LDM coefficients in Eqs.(27) and (34). 

We calculated the various LDM energy coefficients using the analytical expressions 
given in Eq.(28) for the three Skyrme forces SkM*, SkA and S3. In this calculation, only 
the curvature a c and Gauss curvature a Q coefficients have a (3 C dependence. We found 
that both these coefficients are almost invariant with respect to j3 c . Therefore, we have 
used p c = while calculating the energy coefficients, as is normally done. Values obtained 
for the LDM energy coefficients using Eq.(28) are given in Table I for the three Skyrme 
forces. The values obtained here for the various LDM coefficients agree qualitatively well 



with those found in literature ||15|| . In regard to the nuclear curvature, it may be noted 
that the semi-classical ETF models [|18|, |1£|] give a value of about 10 MeV, as against 
the value of about zero, determined |20| from experiments. This is the so-called nuclear 
curvature energy anomaly, which is yet to be resolved. 

Similarly, we calculate the values of compressibility coefficients using Eq.(35). In this 
calculation, we have shown our results for three particular values of f3 c in Table II. The 
ratio | K s /K v | is also tabulated. It can be seen that as f3 c decreases from (3 C = to 
f3 c = —1/2, | K s /K v | also decreses. This behaviour of K s is in agreement with that of 
the pocket model flQfl . Further, we find the scaling model result K s ~ —K v is found to 
be well satisfied for (3 C = —1/3. 

In our study we have also calculated the higher-order coefficients like K c and K for 
different values of (3 C . The value of K c remains more or less unchanged with respect to f3 c . 
Due to this, as f3 c becomes more and more negative, curvature effect becomes as important 
as the surface. This may be due to the fact as f3 c becomes more and more negative, the 
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finite-size effects play a progressively important role in the breathing vibrations of a 
nucleus. Further, from the values of all finite-size compressibility coefficients like K s , K c 
etc given in Table II, it can be seen that (3 C > will lead to a relatively slow converging 
series compared to the other negative values of f3 c . In other words, as f3 c becomes more 
and more positive, | K s /K v | increases. 



It may be mentioned here that in Ref. [13|], we had calculated the coefficients using 
P c = and neglecting terms like e',M _1//3 , e"A~ 2 / 3 etc in the denominator of Eq.(23). 
Now retaining those terms and performing binomial expansion as given in Eq.(24), we 
find that the values of K c and K a obtained with j3 c =0 get substantially modified. For 
reasons discussed in the introduction, we presently study the convergence properties of 
the leptodermous expansion of Ka with f3 c = 0, and also make a comparative study of 
the LDM expansions of energy and compressibility. 

B. Convergence behaviour and Pair effect 

To study the convergence of the LDM expansion of both energy and compressibility, 
it is necessary to evaluate exactly the values of energy E and compressibility Ka corre- 
sponding to symmetric(N=Z) systems, which will then be compared with the sum of all 
the terms given by the right hand side of Eqs.(27) and (34) respectively. We determine 
the total energy numerically using Eqs.(5-10) and then, Ka can be obtained using Eq.(30) 
with j3 c = 0. Here, the diffuseness parameter a is held fixed at its semi-infinite nuclear 
matter(SINM) value. The values so obtained for four representative nuclei in the range 
40 < A < 200 are presented in Tables III and IV for the energy and compressibility re- 
spectively for Skyrme forces SkM*,SkA and S3. It can be seen that in the case of energy, 
both the results agree extremely we\\,i.e. upto about second decimal place, and even, 
in the case of compressibility the agreement is quite good. This result implies that the 
convergence of the LDM expansion of Ka is almost as good as the LDM expansion of 
E/A. 

To understand this inference in a better way, we have presented in Tables V and VI, 
the contributions of the successive terms in the energy and compressibility expansions 
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respectively, for the three forces and, for two representative mass numbers A=40 and 
200. It can be seen that in the case of energy, the successive terms decrease by a factor 
of about 5, giving rise to a rapidly converging series . In contrast to this, we find in 
the case of compressibility, [see Table VI ], the Gauss curvature term KqA" 1 is of the 
same order as the curvature term K C A~ 2 I 3 . In fact, for values of A below about 150 
,KqA~ 1 even overshoots K C A~ 2 / 3 . Same is also found to be true in our calculation for the 
higher-order terms K^A^I 3 and KfA~ 4 ^ 3 . Further, it is interesting to note that K C A~ 2 I 3 
and K$A~ X are almost equal in magnitude, but opposite in sign, which is also the case 
with Kf and K h . Hence, we find that the terms of higher order than the surface one 
in the LDM expansion of Ka cancel pairwise. This pair effect gives rise to a misleading 
conclusion regarding the importance of higher-order terms in the LDM expansion of Ka, 
unless otherwise investigated. Thus, the LDM expansion of Ka in the case of pure bulk 
mode shows an anomalous behaviour, in the sense that KqA" 1 ~ — K C A~ 2 / 3 , in contrast 
to the rapidly converging energy expansion. 

In addition, it may be mentioned here that the convergence properties of the LDM 
expansion of Ka in the case of scaling mode has been widely studied[6-9], going upto order 
of curvature term. In our study, including terms upto Gauss curvature order, we arrive 
at similar conclusions. Although, pair effect is not observed in the scaling case, the rate 
of convergence is still found to be relatively slow as compared to the energy expansion. 

Further, we would like to state that the above results were found to remain valid 
when one uses more realistic ETF functionals including h 4 terms, and with a generalised 
Fermi-function for p{r). In view of this new feature in the convergence behaviour of the 
LDM expansion of Ka, the question arise: whether we can extract the various coefficients 
in the LDM expansion of Ka using a least squares fit analysis, which may manifest the 
effect of correlations amongst the coefficients. 

IV. LEAST SQUARES FIT ANALYSIS 

In the earlier section, we have obtained the exact values of energy E/A and com- 
pressibility Ka numerically and also, the values for the various LDM coefficients in their 
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expansions for SkM*,SkA and S3 forces, and have established their goodness. Now, in 
this section, we would like to determine these coefficients again from a least-squares fit 
to the exact values of E/A and Ka (referred to as synthetic data) for a set of nuclei. A 
comparison of the two sets of values will then demonstrate clearly the reliability of the 
LDM expansions of E/A and Ka for the extraction of nuclear matter properties and, the 
surface properties of finite nuclei. This will also establish the convergence behaviour of 
both the expansions. 

A. For symmetric case 

Here, we attempt to determine all the coefficients in the LDM expansion of Ka per- 
taining to symmteric systems(N=Z) by making a free least squares fit to the numerically 
calculated values of Ka using Eqs. (5-10,30), for 210 nuclei in the mass range 40 < A < 250. 
In Tables VII and VIII, we have presented the results so obtained for E/A and Ka re- 
spectively, for the SkA force only, whose value of K^ lies in between those of SkM* and 
S3 forces. The corresponding error in all the coefficients calculated using the standard 



method |21j are also given. The first row shows the exact values of the coefficients obtained 
in the previous section, which are presented for comparison. 

It can be seen from Tables VII and VIII that, in the case of both energy and com- 
pressibility, the volume term is quite well determined and its value progressively improves 
with the inclusion of higher-order terms. Same is almost true for the surface coeffi- 
cients. And, the higher-order terms such as curvature and Gauss curvature are relatively 
ill-determined. This is presumably due to correlations amongst the various coefficients, 
which is clear from Eqs. (28) and (35). Thus, we find that in the case of symmetric systems 
without the Coulomb effect, the principal coefficients like volume compressibility K v and 
surface compressibility K s can be extracted reliabily by means of a least square fit. 

B. Inclusion of Coulomb and asymmetry effects 

The Coulomb force plays an important role in real nuclei, and also, because of its 
influence in the extraction of K^ from the breathing-mode data as shown by Pearson 
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[0] and Shlomo & Youngblood [ffl, it is worthwhile to consider the effect of Coulomb 
interaction in our present analysis. 

To investigate this, we repeat our calculation including the Coulomb force and asym- 
metry effect, for many realistic nuclei lying within the mass region 40 < A < 250. In this 
calculation of total energy(5), we suppose the neutron and proton density distributions 
at the ground-state to be of the form, 

Pi( T ) = P °1-r,) ! l = n,p (36) 

1 + e a i 

where the various parameters are as defined in Eq.(10), and are determined by energy 
minimisation criteria. Also, for the sake of simplicity, we have used the same value of 
diffuseness parameter for both proton and neutron density distributions, a n = a p . The 
Coulomb energy is given as Ec ou — 0.6Z 2 e 2 / R, where we have considered the nucleus 
to be a uniformly charged sphere of radius R as is normally done|| [8]. One can then 
calculate Ka using the total energy expression(5) including asymmetry and Coulomb 
effects in Eq.(30). 

Now, we make a least squares fit to the so-obtained theoretical values of Ka numbering 
210. The results so obtained from the fit for, the case of compressibility is given in Table 
IX for the SkA force. We find that, in the 4-paramter fit involving K v , K s , Kq u and 
Kp, the extraction of all these coefficients, except for the symmetry compressibility K@, 
are reliable. The discrepancies in K v and K s as compared to the exact values are about 
4% and 12% respectively. With the introduction of a surface-asymmetry K s p term, i.e. 
in the 5-paramter fit, the estimate for Kp improves, however, at the cost of important 
coefficients like K v , K s and Kc ou - To see the effect of a curvature term, we made a 
5-parameter fit (K v , K s , Kcou,Kp & K c ) to the 210 model data on Ka- As it can be 
seen, introduction of a curvature term as a free parameter somewhat spoils the whole fit, 
introducing a maximum of discrepancy of about 40% in K v . 

Therefore, from our analysis using synthetic data, it can be concluded that inspite of 
the correlations among various coefficients, it may be possible to extract K^. However, 
it must be mentioned that even with such synthetic data, inclusion of higher-order terms 
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like K c and K s p as free parameters in the fit, lead to poorer determination of the volume 
coefficient K^, and also other coefficients. Thus, the present theoretical study clearly 
points out the goodness and limitations of the LDM expansion of K A and the extent of 
error inherently present in this approach, in regard to the determination of K^. 

V. DYNAMICAL EFFECTS 

Normally, the convergence properties of the LDM expansion of Ka are mostly studied 
within the scaling model. This is because, only within this model, one can easily relate Ka 
to the experimental data on GMR. However, in general, the mode of density vibrations 
is neither scaling-like nor pure bulk-like. Hence, we need to examine the convergence 
properties of Ka— expansion taking into account this dynamical effect. 

To do so, firstly, we need to relate to the GMR data for any mode of monopole 
vibrations. In other words, we need to find a general empirical relation analogous to 
Eq.(l), and thereby, one can obtain experimental K A values from the GMR data without 
making the scaling assumption. Once K A values are known from the GMR data, one may 
then use Eq.(2) to extract K^. Hence, in the following, we address two aspects: Firstly, 
how can one obtain a general relation between Ka and experimental E gmr ?, Secondly, 
how does Ka behave under the most general conditions? 

In the appendix, we discuss the above mentioned first aspect using the hydrodynamical 
approach. This justifies our study of the Ka— expansion in the generalised situation, i.e. 
without using the scaling assumption. With this, we now focus upon the second aspect, 
i.e. the general behaviour of Ka taking into account the A— dependence of f3 c , in the 
following. 

Within our analytical model, we calculate realistic values of Ka using the general 
expression 

K A (Pc) =K v + K a (p c )A~ 1/3 + K C (P C )A- 2 / 3 + K Cou Z 2 A-^ + K p (3 2 (37) 

where K Cou & Kp are respectively the Coulomb and asymmetry compressibility coeffi- 
cients and the asymmtery paramter (5 = (N — Z) / A. The values for Kcou and K sym 
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obtained with SkM* force are —4.70 MeV and —349.0 MeV. In the calculation of Kcou, 
we have considered only the direct Coulomb term ac = 0.6e 2 /r o . Realistic values of K s 
and K c are obtained as follows. 

For a given nucleus (A,Z), the optimum value of f3 c that will give rise to an excitation 



energy fiu) close to the experimental value can be obtained using[l I 



where (3 is related to A as 



Pc = P - + 1 (38) 



p = 0.685A 1/3 -2.15 (39) 



So, for a given A and the corresponding (3 C obtained from Eq.(37), we have calculated 
K s and K c for the SkM* force using Eq.(38). The values thus obtained for K s and K c 
are plotted as a function of /3 C in Figs.(l) and (2) respectively. It can be seen that as (3 C 
varies from —0.6 to 0.0, K s /K v varies from from about —0.5 to —2. On the other hand, 
variation of K c is not so prominent. Further, it is interesting to note that value of K c 
shows a minimum at about the scaling model value {f3 c = —1/3). An important point 
to be noted is that since /3 C is A-dependent, and K s & K c are [3 C — dependent, the surface 
and curvature compressibility coefficients have a residual mass dependence. To obtain 
the unique A— independent value of K s and K c , one should take the limit A — > oo, 
which leads to (3 C = 0. The asymptotic value of K s thus obtained is nothing but the 
value K s ~ —2K V , as noted earlier [0. And, the true asymptotic value of curvature 
compressibility coefficient is the value of K c obtained with (3 C = plus the contribution 
coming from the f3 c dependence of K s . 

Using the values of K s and K c calculated for the SkM* force with the optimum values 
of j3 c in Eq.(38), we plot in Fig. (3), the so-determined Ka{(3 c ) / K v values as a function 
of A -1 / 3 . For sake of comparison, we have shown in Fig. (4) the values of Ka(Pc)/ Kv 
obtained for three particular values of f3 c . It may be mentioned that, for A > 250, the 
Coulomb force is switched off and N=Z. It can be clearly seen from Fig. (3) that 
shows an ' up - turn ' behaviour as against the nice linear behaviour found for the three 
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particular values of (3 C in Fig. (4). This up-turn behaviour or change in the slope may be 
suggesting the onset of the breakdown of the leptodermous expansion of Ka below mass 
number approximately 120. Indications for such an increasing nature of Ka can also be 
found from the hydrodynamical calculations JIT] . 



In analysing real data, one should indeed expect such an 'up-turn' behaviour as against 
the nice linear behaviour obtained under scaling or pure bulk mode assumption, which 
suggests that higher-order coefficients like K c become important over medium and low 
mass regions. Because of this 'up-turn' behaviour , it may be difficult to consistently 
determine all the parameters in Eq.(2) from a fit to the presently available few tens of 
data on GMR, which in-turn shall impair the extraction of the important quantity, K^. 

VI. CONCLUSION 

In conclusion, the LDM expansion of Ka in the pure bulk mode, shows an anomalous 
convergence behaviour due to pair effect, as compared to the widely studied scaling mode. 
It is also found that the Ka— expansion in both the cases is relatively quite slow, as 
compared to the energy expansion. However, K^ can be reliably extracted in these 
specific cases. In realistic situations, one also encounters modes of density vibrations, 
other than the scaling and pure bulk modes, depending upon the mass region under 
consideration. When this dynamical effect is taken into account, the nuclear compression 
modulus shows an 'up-turn' behaviour below mass number about 120, suggesting the 
inapplicability of the LDM expansion of Ka over this mass region. 
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APPENDIX 

Here an attempt is made to obtain a generalised empirical relation between the nuclear 
compression modulus Ka and the experimental breathing- mode energies E gmr . Until now, 



a well-defined relation E gmr = J % Ka/ (m < r 2 >), is available only in the case of scaling 
model. 

Brack and Stocker[|TT|1 using a variational hydrodynamical approach found that, the 
mass paramter B(j3 c , A) shows a regular behaviour with respect to both coupling paramter 
f3 c and mass number A. Using their results displayed in Table 1 of Ref . Jll[ , we have plotted 



B((3 C , A)/m < r 2 > as a function of A~ 1 ^ in Fig. (5). It can be seen that the general mass 
paramter B(/3 C ,A) expressed in terms of the scaling model value m < r 2 > varies quite 
smoothly with respect to A. Hence, B(j3 c ,A) can be expressed as 

B(P C ,A) ~ m<r 2 >f(A) 

where f(A) = Co + C\A~ 1 ^ + • ■ ■ can be a polynomial in A^ 1 ^ 3 . The crucial function f(A) 
can be determined from hydrodynamical calculations, which will be reported elsewhere. 
Once f(A) is known, one can determine Ka from GMR data using the relation 



K A (f3 c 



K A ((3c 



\B((3 C ,A) \\m<r 2 >-f(A) 

Thus, a generalised expression relating Ka and E gmr , analogous to the well-known scaling 
relation, is proposed. This proposition is on sound footing as the hydrodynamical studies 
are successful in describing the GMR data. 
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TABLE CAPTIONS 



Table I: Values of the various coefficients in the LDM expansion of the energy(27) for 
the Skyrme forces SkM*, SkA and S3 using the local ETF functional with h 2 terms and 
a pure Fermi-function for the density distribution(3). All quantities are in MeV. 

Table II: Values of the various coefficients in the LDM expansion of the compressibil- 
ity (34) for the Skyrme forces SkM*, SkA and S3 using the local ETF functional with h 2 
terms and a pure Fermi-function for the density distribution . Three values of f3 c have 
been used. All quantities are in MeV. 

Table III: Comparison of the exact values of the total energy per nucleon E/A (5) with 
the analytically determined values using the LDM expansion of E/A, given by Eq.(27), 
for four representative mass numbers A, and for the Skyrme forces SkM*, SkA and S3. 
All quantities are in MeV. 

Table IV: Comparison of the exact values of the finite nuclear compression modulus 
Ka (30) obtained using f3 c = with the analytically determined values using the LDM 
expansion of Ka , given by Eq.(34), for four representative mass numbers A, and for the 
Skyrme forces SkM*, SkA and S3. All quantities are in MeV. 

Table V: Values of the different terms contributing to the total energy per nucleon 
obtained using the Skyrme forces SkM*, SkA and S3, for two representative mass numbers 
A. All quantities are in MeV. 

Table VI: Values of the different terms contributing to the nuclear incompressibility 
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obtained using the Skyrme forces SkM*, SkA and S3 with (3 C = 0, for two representative 
mass numbers A. All quantities are in MeV. 

Table VII: Values of the parameters obtained from a least-squares fit to the exact values 
of energy per nucleon E/A obtained for symmetric systems using the SkA force . The 
number of data points used is 210, in the mass region 40 < A < 250 . The first row 
gives the exact values for the various coefficients obtained in our analytical model. All 
quantities are in MeV. 

Table VIII: Values of the parameters obtained from a least-squares fit to the exact values 
of nuclear incompressibility Ka obtained for symmetric systems using the SkA force with 
(3 C = 0. The number of data points used is 210, in the mass region 40 < A < 250 . 
The first row gives the exact values for the various coefficients obtained in our analytical 
model. All quantities are in MeV. 

Table XI: Values of the parameters obtained from a least-squares fit to the exact values of 
nuclear incompressibility Ka obtained for asymmetric systems with Coulomb interaction, 
using the SkA force and taking f3 c = 0. The number of data points used is 210, in the mass 
region 40 < A < 250 . The first row gives the exact values for the various coefficients. 
Value of K s p obtained within the scaling model is taken from Ref. ||. All quantities are 
in MeV. 
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FIGURE CAPTIONS 
FIG. 1 Values of the mass paramter B(f3 C) A) expressed in terms of m < r 2 >, obtained 
in a hydrodynamical calculation with SkM* force is plotted versus A -1 / 3 , where m is 
the nucleon mass and < r 2 > is the root mean square radius. 

FIG. 2 Values of the ratio of the surface compressibility coefficient K s to the nuclear 
matter incompressibility K v obtained in our analytical model using SkM* force is shown 
as a function of the coupling paramter f3 c . 

FIG. 3 Same as Fig. 2, but for the curvature compressibility coefficient K c . 
FIG. 4 Values of the ratio of finite nuclear compression modulus to K v obtained 
including the dynamical effect (A— dependence of /3 C ) is shown versus A~ 1/3 . The force 
used is SkM*. 

FIG. 5 Values of the ratio of finite nuclear compression modulus Ka to K v obtained for 
three particular values of f3 c is shown versus A^ 1 ^ 3 . The force used is SkM*. 
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Table I 





SkM* 


SkA 


S3 


a v 


-15.79 


-16.01 


-15.87 


a s 


19.08 


19.95 


18.90 


a c 


10.24 


9.67 


6.87 


a 


-12.21 


-11.42 


-7.24 



Table II 



Force 


Pc 


K s (MeV) 


K c (MeV) 


K (MeV) 


K s /K v 


SkM* 





-406.1 


-109.9 


568.1 


-1.87 


(K v = 217 MeV) 


-1/3 


-231.0 


-129.3 


138.1 


-1.07 




-1/2 


-129.1 


-118.0 


-69.1 


-0.6 


SkA 





-484.6 


-123.7 


595.0 


-1.84 


(K v = 263 MeV) 


-1/3 


-295.8 


-145.7 


167.8 


-1.12 




-1/2 


-186.4 


-138.0 


-43.8 


-0.71 


S3 





-570.3 


-114.1 


452.5 


-1.60 


(K v = 356 MeV) 


-1/3 


-389.6 


-140.1 


156.2 


-1.10 




-1/2 


-285.1 


-142.7 


5.26 


-0.8 
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Table III 



A 


SkM* 


SkA 


S3 




Exact Eq.(27) 


Exact Eq.(27) 


Exact Eq.(27) 


40 


-9.72 -9.64 


-9.70 -9.64 


-9.96 -9.94 


100 


-11.35 -11.33 


-11.39 -11.38 


-11.56 -11.55 


150 


-11.93 -11.92 


-12.00 -12.00 


-12.12 -12.12 


200 


-12.30 -12.39 


-12.38 -12.37 


-12.47 -12.47 



Table IV 



A 


SkM* 


SkA 


S3 




Exact 


Eq.(34) 


Exact Eq.(34) 


Exact 


Eq.(34) 


40 


103.4 


102.7 


126.9 125.9 


191.3 


190.3 


100 


130.2 


129.7 


159.5 159.1 


232.3 


231.9 


150 


140.4 


140.1 


172.0 171.7 


247.4 


247.1 


200 


147.1 


146.8 


180.0 179.8 


257.1 


256.9 



27 



Table V 



Force 


A 


a v 


asA- 1 ' 3 


a c A- 2 ^ 


aoA- 1 




40 


-15.79 


5.58 


0.8755 


-0.3052 


SkM* 


200 


-15.79 


3.26 


0.2994 


-0.0611 




40 


-16.01 


5.83 


0.8268 


-0.2855 


SkA 


200 


-16.01 


3.41 


0.2828 


-0.0571 




40 


-15.87 


5.53 


0.5874 


-0.1810 


S3 


200 


-15.87 


3.23 


0.2009 


-0.0362 



Table VI 



Force 


A 


K v 




K C A- 2 / 3 


KoA- 1 




40 


216.6 


-118.7 


-9.40 


14.20 


SkM* 


200 


216.6 


-69.4 


-3.21 


2.84 




40 


263.3 


-141.7 


-10.58 


14.88 


SkA 


200 


263.3 


-82.9 


-3.62 


2.98 




40 


355.5 


-166.8 


-9.76 


11.31 


S3 


200 


355.5 


-97.5 


-3.34 


2.26 
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Table VII 



No. of 
Para. 


a v 


a s 


a c 


a 




-16.01 


19.95 


9.67 


-11.42 


2 


-16.2±0.0003 


22.1±0.002 






3 


-16.2±0.002 


22.2±0.02 


-0.30±0.04 




4 


-16.0±0.0001 


19.7±0.001 


11.6±0.007 


-18.2 ± 0.01 



Table VIII 



No. of 
Para. 


K v 


K s 


K c 


K 




263.3 


-484.6 


-125.8 


604.4 


2 


256.4±0.14 


-447.7±0.69 






3 


269.0±0.07 


-569.1±0.65 


282.1±1.5 




4 


263.7±0.24 


-494.0±3.4 


-67.1±15.8 


531.2±24.0 
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Table IX 



K v 


K s 




Kcoul 


K c 


K sf s 


263.3 
252.3±2.0 
237.5±5.8 
368.7±17.3 


-484.6 
-428.4±5.6 
-391.2±14.7 
-1108.4±100.5 


-441.1 
-240.8±7.0 
-394.6±56.8 
-423.3±28.0 


-5.14 
-4.7±0.19 
-3.28±0.55 
-9.4±0.71 


-125.8 
1126.5±166.3 


~ 875 
1009.8±370.0 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9505031vl 



